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ABSTRACT 

We fit the broad-band X/7-ray spectrum (0.8-1000 keV) of the accreting black hole XTE 
J1550-564 in the very high state. The quasi-simultaneous data from A5CA, RXTE and OSSE 
show that the disc is strongly Comptonized, with a high energy tail extending out to several 
hundreds keV. However, inverse Compton scattering by a purely thermal or purely power-law 
electron distribution cannot explain the observed spectrum. Instead the data require a hybrid 
distribution, with both thermal and non-thermal electrons scattering the disc photons. This is 
very similar to the electron distribution inferred for other high and very high state black hole 
binaries, showing that it is a generic feature of high mass accretion rate black holes. 

Key words: accretion, accretion discs - X-rays: individual: XTE J1550-564- X-rays: bina- 
ries - radiation mechanisms: non-thermal 



1 INTRODUCTION 

Galactic black hole candidates show variety of spectra at X-ray 
and 7-ray energies. These can be roughly described in terms of 
a soft, quasi-thermal component from the accretion disc together 
with a power-law-like tail extending to much higher energies. Five 
'canonical' spectral states have been identified historically, based 
on the luminosity and spectral shape (see e.g. Esin, McClintock & 
Narayan 1997 and references therein). The very high state is gen- 
erally seen at the highest luminosities, where the soft thermal com- 
ponent (with temperature of 0.5-1 keV) and power-law tail (photon 
spectral index T = 2-3) have comparable power. In the high state 
the hard tail is much weaker or sometimes even not present. The 
intermediate state is rather like the very high state, but at lower 
disc temperature and luminosity (Belloni et al. 1996; Mendez & 
van der Klis 1997), while the low state is somewhat different, be- 
ing characterised by a hard (F < 2) power-law-like spectrum with 
a very weak, low temperature disc component (Bahicinska-Church 
et al. 1995; di Salvo et al. 2001). At very low luminosities, transient 
sources go into the quiescent state, which is a very low luminosity 
version of the low state (Kong et al. 2000; Sutaria et al. 2002). The 
high and very high states are often described as soft states, while the 
low state is commonly referred to as a hard state. Persistent sources 
like Cyg X-1 and LMC X-1 spend most of their time in one state, 
however transients like Nova Muscae can progress through a se- 
quence of all five states during an outburst (Ebisawa et al. 1994). 
See Done & Gierlinski (2003) for an overview of these spectral 
states. 

The power-law-like tails extending to X/7-ray energies are 
commonly explained by Compton up-scattering of the soft seed 
photons from the accretion disc by a population of high-energy 
electrons. In the low/hard state the observed hard spectra roll over at 
100-200 keV, indicating that the electrons are mostly thermal (with 



typical temperatures of '-^ 100 keV) and have optical depth of or- 
der unity (e.g. Gierlinski et al. 1997; Zdziarski et al. 1998; Frontera 
et al. 2001). By contrast, in the soft states the steeper tail extends 
unbroken to MeV energies (Grove et al. 1998), indicating that the 
electron energy distribution is predominantly non-thermal. More 
detailed spectral studies of the very high- and high-state spectra 
show that the tail is best described by Comptonization from a com- 
plex electron distribution, where there is both a low temperature 
(~ 10 keV) thermal electron distribution (Zycki, Done & Smith 
2001; Kubota, Makishima & Ebisawa 2001) together with non- 
thermal power-law electrons (Gierliiiski et al. 1999; Zdziarski et 
al. 2002). A two-component electron distribution could arise fairly 
naturally from non-thermal electron acceleration regions powered 
by magnetic reconnections above a disc. Low-energy electrons cool 
preferentially by Coulomb collisions leading to a thermal distri- 
bution while the high-energy electrons cool by Compton scatter- 
ing, preserving a non- thermal distribution (Coppi 1999). Alterna- 
tively, the thermal and non-thermal electrons could be spatially dis- 
tinct, e.g. magnetic reconnection above the disc can produce a non- 
thermal electron distribution, while overheating of the inner disc 
produces the thermal Comptonization (Kubota et al. 2001). 

The soft and hard states are clearly distinct in their spectral 
properties (e.g. Zdziarski et al. 2002). Transition between the states 
is compatible with being driven by a changing accretion geome- 
try, where the low state has an accretion disc only at large radii, 
with the inner disc being replaced by a hot, optically thin, geo- 
metrically thick flow (e.g. Poutanen, Krolik & Ryde 1997). With 
increasing accretion rate the disc moves inwards increasing its tem- 
perature and the number of seed photons cooling the inner flow so 
the Comptonized spectrum becomes softer. Eventually, when the 
inner flow becomes optically thick it collapses, the inner disc re- 
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places the hot flow and the source reaches one of the soft states 
(see e.g. Esin et al. 1997). 

This review of spectral properties draws heavily from the very 
few broad bandpass X/7-ray observations of black holes. To con- 
strain both the disc and hard component requires simultaneous data 
over a broad energy range, ideally ~ 0. 1-1000 keV. While both the 
Phoswich Detector System on board BeppoSAX and High Energy 
X-ray Timing Experiment (HEXTE) on board Rossi X-ray Timing 
Explorer (RXTE) can give data out to ~ 100-200 keV, the Ori- 
ented Scintillation Spectrometer Experiment (OSSE) on the Comp- 
ton Gamma Ray Observatory (CGRO) is unique in being able to 
extend the bandpass to higher energies. This is especially impor- 
tant in soft states, as there the hard component is known to ex- 
tend to the highest energies. However, there are only very few 
broad-band observations of black hole candidates in these spectral 
states: GRS 1915+105 (Zdziarski et al. 2001; Ueda et al. 2002), 
GRO J1655-40 (Tomsick et al. 1999), Cyg X-1 (Gierlinski et 
al. 1999; McConnell et al. 2002). However, both GRO J1655-40 
and GRS 1915+105 are thought to be extreme Kerr black holes due 
to their high accretion disc temperatures (e.g. Zhang, Cui &. Chen 
1997), so there are no very high state broad-band spectra of 'nor- 
mal' black holes. 

In this paper we remedy the situation by analysing the broad- 
band simultaneous RXTE/OSSE data of the X-ray black hole tran- 
sient XTE J 1550-564, and derive constraints on the electron dis- 
tribution and accretion flow geometry. We show that the hard tail 
is remarkably similar to that seen in the extreme Kerr, microquasar 
black holes, so it is most probably a generic feature of high mass 
accretion rate flow onto a black hole rather than being associated 
with spin or radio jets. 



2 THE SOURCE 

The X-ray nova XTE J 1550-564 was discovered during its spectac- 
ular outburst in 1998 (Smith et al. 1998). Two weeks into the out- 
burst it reached a peak flux of 6.8 crab (1.6 x 10~^ erg s~^ cm~^) 
in 2-10 keV energy range, making it the brightest X-ray transient 
observed so far by RXTE. The X-ray (2-10 keV) lightcurve from 
All-Sky Monitor (ASM), is more complex than a fast rise, expo- 
nential decay, having a bright flare followed by two plateaus sep- 
arated by a weak minimum. Such lightcurves are not uncommon 
(see e.g. Chen, Shrader & Livio 1997), and are probably connected 
to enhanced mass transfer from the irradiated companion star (Esin, 
Lasota & Hynes 2000). The complex X-ray spectral evolution dur- 
ing the outburst is shown in Fig. Q by comparing the 2-10 keV 
lightcurve from ASM to the harder X-ray (20-100 keV) lightcurve 
from CGRO Burst And Transient Source Experiment (BATSE). 
More details of this spectral evolution are given in Sobczak et 
al. (2000), Homan et al. (2001), Wilson & Done (2001) and Kubota 
et al. (2003). 

Observations of the optical companion (Orosz et al. 2002) led 
to an estimate of the compact object mass, 8.4-10.8 M0, confirm- 
ing its black hole nature. The inclination angle of the binary is quite 
high, 67°-77.4° (both estimates are given with la errors, the 3a 
ranges are 6.8-15.6 Mq and 55°-79°, respectively). The distance 
is less well constrained by their data on the optical companion, with 
D = 3.0-7.6 kpc (best fit value of 5.3 kpc). Throughout this paper 
we assume M = WMq, i = 70° and D = 5.3 kpc. 
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Figure 1. BATSE and ASM lightcurves of 1998 outburst of XTE J1550- 
564. Vertical line marks the OSSE viewing period 729.5. analysed in this 
paper together with simultaneous RXTE observations. The ASCA observa- 
tions was made one day before the OSSE viewing period 729.5 had started. 



3 THE DATA 

In this paper we analyse the contemporaneous X/7-ray spectra of 
XTE J1550-564 from ASCA, RXTE and OSSE in a period from 
1998 September 23 to 1998 October 6 (MJD 51080-51092) i.e. 
in the first part of the outburst, in the plateau just after the strong 
flare. There are RXTE observations over this entire period, while 
the ASCA data are limited to a snapshot at the start (September 
23/24) which does not quite overlap with the OSSE data (between 
September 25 and October 6). The times of our data are shown on 
the soft and hard X-ray lightcurves of the outburst in Fig.Qwhile a 
detailed observation log is presented in TableQ 

We used the standard product ASCA Gas Imaging Spectrome- 
ter (GIS) spectrum from detector 2, with response matrix version 
4.0, corrected for dead time (Makishima et al. 1996). Since the 
source is very bright we fit the GIS2 spectrum alone to avoid possi- 
ble cross-calibration uncertainties between GIS2 and GIS3. There 
is a simultaneous RXTE pointing (observation number 30191-01- 
10-00), but we treat the ASCA data separately to avoid the signifi- 
cant cross-calibration issues between RXTE and ASCA (e.g. Done, 
Madejski, Zycki 2000). 

The OSSE data for viewing period 729.5 extends between 
September 25 and October 6 (Fig.0. We use a high-level product 
spectrum extracted for detectors 1, 2, 3 and 4, for the whole viewing 
period. Due to rather large OSSE field of view (1 1° x 3.8° ; Johnson 
et al. 1993) there is always a danger of a background source con- 
taminating the spectrum. Therefore, we compare the OSSE spec- 
trum with the simultaneous HEXTE spectrum. Since HEXTE has 
much smaller field of view (~ 1°; Rothschild et al. 1998), contami- 
nation by a background source is much less probable. The excellent 
agreement between the HEXTE and OSSE spectra (see Section 5.2) 
make contamination by a background source rather unlikely. 

RXTE observed the source 15 times during OSSE viewing pe- 
riod 729.5, with observation archival numbers (obsids) between 
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Table 1. Log of obsei-vations. 



Instrument 


Detectors 


Observation ID 


Time (MJD) 


Live time (s) 


Count rate (cnt s ^) 


ASCA CIS 


2 


15606010 


51079.964-51080.500 


25156 


93.30±0.06 


RXTE PCA 


01, layer 1 


30191-01-12-00 ... -23-00 


51082.003-51091.790 


45664 


5891±6 


RXTE HEXTE 





30191-01-12-00 ... -23-00 


51082.003-51091.790 


14442 


287.8±0.2 


CGRO OSSE 


1234 


729.5 


51081.640-51092.561 


496313 


12.2±0.2 
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Figure 2. RXTE/PCA colour-colour diagram of the 1998 outburst of XTE 
J 1550-564. The colours were corrected for the PCA gain changes (by the 
same method as in Gierlitiski & Done 2002) and rebinned for clarity. Data 
points in black denote the RXTE observations simultaneous with OSSE 
viewing period 729.5. 



30191-01-12-00 and 30191-01-23-00. To assess the source X-ray 
variability over this 10 day period we create a PCA colour-colour 
diagram, using a method described by Gierlitiski & Done (2002). 
The diagram, shown in Fig.|2| contains all the PCA data collected 
during the 1998/1999 outburst, and observations simultaneous with 
OSSE are marked in black. There is a shift in colours during this pe- 
riod, though it is small compared to the spectral variation seen dur- 
ing the whole outburst. This particular part of the colour-colour dia- 
gram corresponds to the very high spectral state (Done & Gierliiiski 
2003). 

We extract the PCA and HEXTE spectra using FTOOLS 5.2 for 
each of the aforementioned obsids, and also create a single spec- 
trum for each instrument averaged over the whole OSSE viewing 
period. For the PCA, we extract the spectra from detectors and 1, 
top layer only (see Wilson & Done 2001), correcting for the photon 
pile-up (Tomsick et al. 1999). For HEXTE we extract spectra from 
cluster only. 



4 SPECTRAL MODELS 

As outlined in the introduction, the X-ray spectrum is probably 
formed by Compton scattering of the accretion disc seed photons. 
While an exponentially cutoff power law is often used as an ap- 
proximation to a thermally Comptonized spectrum, this becomes 
very inaccurate if the observed bandpass approaches the typical 
seed photon or electron energies (e.g. Done, Zycki & Smith 2002; 



Zdziarski et al. 1998). Since we have an extremely broad bandpass 
we use instead the best available Comptonization models, which 
properly calculate the shape of the high-energy rollover for ther- 
mal electrons, and include the low-energy cutoff close to the seed 
photon energies. There are currently two of these Comptonization 
codes available for spectral fitting, COMPPS (Poutanen & Svensson 
1996) and EQPAIR (Coppi 1999). Both of these can calculate the 
spectra from thermal as well as non-thermal Comptonization. 

The COMPPS code computes Comptonization in various spe- 
cific geometries (we assume spherical) from a given electron distri- 
bution of total optical depth r. This is split between a thermal dis- 
tribution (Maxwellian, with temperature Te) which is truncated at 
some Lorentz factor 'ymin, above which it has a non-thermal [power 
law, A''(7) oc 7^^] distribution up to Lorentz factor 7max- 

The EQPAIR code is rather different in that it computes the 
electron distribution self-consistently, assuming that the electrons 
are injected with a power-law distribution of index Pinj. The rate 
of Compton cooling is larger at higher energies, so the equilibrium 
distribution is softer than the injected one with index p = Finj + 1 
(e.g. Blumenthal & Gould 1970). The electrons can also thermal- 
ize through Coulomb collisions, producing a hybrid, thermal/non- 
thermal distribution. However, the code also allows some frac- 
tion of the total power to be used to heat the thermal distribution 
directly. The microscopic processes are calculated in the plasma 
cloud of assumed optical depth of ionization electrons, rion. The 
total optical depth, r — rion + Tpair, is calculated self-consistently, 
including the additional optical depth in electron-positron pairs, 
Tpair, produccd by photon-photon collisions. The plasma proper- 
ties depend on its compactness 

C 



7^' 



(1) 



where £ is the power supplied to the plasma, TZ is its characteristic 
size, (Tt is the Thompson cross section and me is the electron mass. 
The power supplied to the electrons is parameterized as a hard com- 
pactness, while the power in soft seed photons is given as a soft 
compactness, £s. The hard compactness is split into non-thermal 
electron acceleration and direct thermal heating, denoted Inth and 
£th, respectively, where = inth+ ^th- The electron temperature, 
Te, is calculated self-consistently from the balance between cooling 
and heating. We note that even when the whole power is supplied 
to non-thermal electrons (i.e. £th = 0), the Coulomb cooling can 
dominate at lower energies over Compton cooling, so the electrons 
will form a hybrid distribution. The model parameters are £s,£h/£s, 
/4, Finj and 

' ion • 

Both models include the reflected continuum computed from 
the incident continuum using the Green's functions of Magdziarz 
& Zdziarski (1995). This is relativistically smeared (Fabian et 
al. 1989), but does not include the self-consistently produced iron 
Ka line, so we include this separately using the DISKLINE model 
(Fabian et al. 1989), with relativistic smearing set equal to that of 
the reflected continuum. These are parameterized by the solid an- 
gle of the reflecting medium, Q, its ionization, ^, and the inner disc 
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radius, i?rcf- The rest-frame energy of the line is denoted as Eunc 
and its equivalent width as EW. 

For purely thermal Comptonization we also use a more 
approximate (but faster) THCOMP code (Zdziarski, Johnson & 
Magdziarz 1996), which is based on Kompaneets (1956) equation. 
This assumes that the photons diffuse in both time and space (Sun- 
yaev & Titarchuk 1980), which becomes inaccurate for tempera- 
tures above ~ 50 keV, as such photons lose a large fraction of their 
energy in Compton down-scattering. 

In all of these models we assume that the seed photons 
for Comptonization have a multi-colour disc blackbody (model 
DISKBB; Mitsuda et al. 1984) distribution, with temperature 
^sccd — iskbb. We allow for additional disc photons which do 
not pass through the Comptonizing region, e.g. corresponding to a 
fraction of the disc not covered by a patchy corona. For calculating 
inner disc radius we take into account all disc photons, both seed 
and non-seed. In figures we plot only the disc photons which reach 
the observer unscattered. 



5 RESULTS 

For spectral fitting we use the X-ray spectral fitting package XSPEC 
version 11 (Amaud 1996), with our own implementation of the 
Comptonization models described above. The error of each model 
parameter is given for a 90 per cent confidence interval (Ax^ = 
2.7). We use the 0.9-10 keV GIS data, 3-20 keV PCA data, 20- 
150 keV HEXTE data and 50-1000 keV OSSE data. The relative 
normalizations of the PCA, HEXTE and OSSE are uncertain, so 
we allow them to be free parameters in simultaneous spectral fits. 
We base fluxes on the PCA normalization. 

5.1 ASCA 

The near-simultaneous ASCA observation gives us an excellent op- 
portunity to study the soft X-ray part of the XTE J1550-564 spec- 
trum, not covered by RXTE and OSSE. From the PCA data we 
know that the X-ray spectrum did not change much between ASCA 
and OSSE observations, so we can apply the seed photon temper- 
ature, Tsccd, and Galactic absorption, Nh, obtained from A5CA to 
the RXTE/OSSE spectral fits. 

We use all three Comptonization models to check that the re- 
sults do not depend on the details of the models used. The ASCA 
data are not sensitive to the details of the high energy spectrum, so 
we assume that the Comptonizing electrons are thermal with elec- 
tron temperature fixed at kT^ = 30 keV (COMPPS and THCOMP) or 
equivalently rion = 2 for EQPAIR. 

Our model consists of this thermal Comptonization (together 
with its reflected continuum and associated iron line emission), 
and a multicolour disc blackbody. All three models for the Comp- 
tonized continuum give almost identical values of Nh and Tsccd 
(Table|2j. We check that our choice of Te (or equivalently rion) did 
not affect these values, and even replacing the disc spectrum by a 
single-temperature blackbody gives only a slightly lower column 
of about 6.0 X 10^^ cm~^ 

This result is not very different from previous findings. Jain et 
al. (1999) found Nh ~ 9 x 10^^ cm"^ from optical reddening in 
the direction of XTE J 1550-564. Kong et al. (2002) fit the Chan- 
dra spectra of XTE J 1550-564 in quiescence with various models 
and found Nh between 3 (blackbody) and 9 x 10^^ cm^^ (power 
law). Thought it not clear what spectral shape should we expect 
from a black hole binary in a quiescence, Tomsick, Corbel & Kaaret 



Table 2. Fits of a multicolour disc blackbody plus thermal Comptonization 
(and its associated reflection spectrum and iron line) to the ASCA spectrum. 
The seed photons for Comptonization are assumed to be from the disc, with 

^sccd — Tin- 



Model 


Nh (10^1 cm-2) 


fcTdiskbb (keV) 




THCOMP 
COMPPS 
EQPAIR 


f. 40+0.09 

6 49+°-0^ 
"•^''-0.12 

"■^■^-0.10 


0.5291°°- 
0.528«°1? 
0.531«:°- 


209.0/184 
206.7/184 
207.4/184 



(2001) point out that the blackbody model leads to rather unphysi- 
cal size of the source, ~ 0. 1 km. Raymond-Smith, bremsstrahlung 
and power law models with Chandra data all give Nh consistent 
with our result. We note that a much larger column of 2 x 10^^ 
cm~^ found by Sobczak et al. (2000) from the PCA fits is due to 
lack of sensitivity of the PCA below ~ 3 keV. And indeed, when 
we allow Nh to be a free parameter in our RXTE/OSSE fit (see 
below, Section|53}, it yields Nh ~ 2.5 x 10^^ cm~^. 

Strongly ionized reflection is significantly detected in all these 
models (see also Wilson & Done 2001), but rather surprisingly the 
amount of relativistic smearing is small. From all three models we 
find very high inner disc radius, typically Ricf > 200-Rg (where 
Rg = GM/c^ is the gravitational radius). However, closer investi- 
gation of the residuals shows that this is driven by narrow features 
in the spectrum around the iron line energy. A disc which extends 
much further into the gravitational potential is consistent with the 
data if there is a narrow absorption line at ~ 7.2 keV with 40 eV 
equivalent width. Such narrow features are seen from other X-ray 
binary systems viewed at high inclination, as GRO J 1655-40 (Ueda 
et al. 1998) or GRS 1915+105 (Kotani et al. 2000), where it is inter- 
preted as resonance absorption from iron, but the line energy we ob- 
tain in XTE J1550-564 is slightly too high. A P Cygni line profile, 
as seen in Cir X-1 (Brandt & Schulz 2000), gives an equally good 
fit to the XTE J 1550-564 data, but with emission and absorption 
line energies ('-^ 6.9 and 7.0 keV, respectively) which are physi- 
cally consistent with the resonant transition in Fe XXVI. However, 
the seed photon energy and absorption column are unchanged by 
these different models of the iron line structure. Thus the derived 
Nh and Tsccd are fairly robust. 

5.2 OSSE and HEXTE 

The OSSE spectrum extracted from the entire viewing period give 
an excellent fit (x^ jv = 7.5/19) to a simple power law with pho- 
ton spectral index F — 3.09 ± 0.04. There is no evidence for any 
high-energy cutoff in the data with a lower limit on an exponential 
rollover energy of 990 keV. 

The HEXTE data are evenly spread throughout the OSSE ob- 
servation, so we can use the better signal to noise of the HEXTE 
instrument to look for spectral variability. The PCA data clearly 
shows that at lower energies there is indeed a small spectral change 
throughout this time, with root-mean-square variation in the soft 
and hard PCA colours of 2.9 and 2.6 per cent, respectively (Fig.|2}. 
We fit the HEXTE spectrum from each RXTE obsid by a cutoff 
power law and Fig.|3|shows the resulting values of the e-folding en- 
ergy, spectral index and flux. The spectral shape is consistent with 
remaining constant (with F ~ 2.4 and Ec ~ 80 keV) despite the 
20-150 keV flux decreasing by ~ 15 per cent. Hence we are justi- 
fied in co-adding all the HEXTE and OSSE data accumulated over 
this period to form a single spectrum from each instrument. 
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Figure 3. Variability of XTE J1550-564 in 20-150 keV range during OSSE 
observation. The HEXTE spectra were fitted by the cutoff power law and 
the panels show the e-folding energy, E^, photon spectral index, F. and 20- 
150 keV model flux in 10~^ erg cm~^ units. Horizontal lines show 
average values of these parameters. The last observation for which only 
lower limit on Ec exists has been excluded from average of Ec- 
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Figure 4. HEXTE and OSSE residuals. The data were fitted by a thermal 
Comptonization plus power law model. The relative normalization between 
HEXTE and OSSE was free. 

Plainly, the power law and exponential cutoff models have 
very different parameters in the HEXTE and OSSE spectra, so a 
joint fit with this model gives a very poor x^- Fig.|4|shows the resid- 
uals of a joint fit of the two spectra to a phenomenological model of 
a power law plus thermal Comptonization model. The fit is reason- 
ably good (x^/u = 65.8/54), but more importantly, Fig. |4| clearly 
shows that the HEXTE and OSSE data are in very good agreement 
in their region of overlap. There are no big cross-calibration prob- 
lems, or problems with a contaminating source in the large OSSE 
field of view. 



5.3 RXTE and OSSE 

Finally, we fit the averaged PCA, HEXTE and OSSE data together 
with the proper Comptonization models. We fix Nh ~ 6.5 x 10^^ 
cm~^ and fcTscod = fcTdiskbb = 0.53 keV from the ASCA results 
(Section lsHl . We estimate the inner radius from the DISKBB model, 
Jidiskbb, assuming distance of D = 5.3 kpc, inclination angle of 
i = 70°. We correct it for spectral hardening due to scattering, as- 
suming the ratio of colour to effective temperature of 1.8 (Shimura 
& Takahara 1995; Merloni, Fabian & Ross 2000). We also cor- 
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rect it for the inner torque-free boundary condition (Gierlihski et 
al. 1999). 



5.3. 1 Maxwellian and power-law electrons 

Purely thermal Comptonization (a model referred to as TH) cannot 
fit the observed spectrum. A model consisting of the disc black- 
body, thermal Comptonization (any code described in Sec.|4} and 
its reflection/line gives /v ^ 5. The best fit is shown in Fig.lsjj 
and clearly shows that the spectral curvature in the lower energy 
data, which requires a fairly low temperature (fcTe ~ 80 keV, r ~ 
0.5) thermal Compton component, is strongly inconsistent with the 
higher energy data. In principle, a very high electron temperature 

100 keV) could fit the high-energy tail, but to have this to- 
gether with a fairly steep overall spectrum requires a very low op- 
tical depth of the Comptonizing medium (r < 0.1). The spectrum 
from Comptonization in such an optically thin plasma consists of 
distinct humps from subsequent scattering orders, which are not 
present in the data. Moreover, the fraction of unscattered disc pho- 
tons would be much larger than seen in the data. 

The high-energy tail in the data can be reproduced in the 
model by adding a power law. This gives much better, though still 
statistically unacceptable fit ix^ /i^ = 161/93). This is because the 
high-energy data requires fairly hard power law (F ~ 2.3 from the 
fit) which extends down to lower energies and strongly contributes 
to the spectrum at soft X-rays, distorting it considerably. 

Interestingly, a purely non-thermal, power-law electron distri- 
bution (model PL) cannot fit the data either (x^ /f ~ 4; Fig.!?!*)- 
This is because the spectrum from multiple Compton scatterings 
on power-law electrons is not a power law. The steep spectrum 
requires a soft electron index (best fit value of p = 6.3 from 
COMPPS). This means that the average electron energy is low, so 
photons must be scattered many times before they reach high en- 
ergies (> 100 keV). In such multiple scatterings they reach the 
Klein-Nishina regime where the cross-section for the scattering de- 
creases. This creates a break in the spectrum at a few hundred keV 
(Ghisellini 1989), so underpredicting the higher energy data above 
~ 300 keV. Power-law electrons with very small optical depth 
would be able to produce an unbroken power-law photon spectrum 
on a single scattering. This, however, would give too large fraction 
of unscattered disc photons again. 

Plainly, the data require a more complex continuum. Now we 
use a model PLTH, consisting of two Comptonizing regions (both 
computed with COMPPS): one with a Maxwellian (thermal) elec- 
tron distribution, the other with power-law electrons. Both these 
Comptonize the seed photons from the disc independently. We ob- 
tain a very good fit with x^/'^ = 65.9/94 (see Table|3j- The spec- 
trum at lower energies is dominated by the thermal component with 
temperature about 30 keV, while the high-energy tail is created 
by up-scattering seed photons on optically thin (r ~ 0.3), steep 
power-law (p ~ 4.5) electrons, where the break from the Klein- 
Nishina is masked by the thermal Compton component (Fig.|6|j). 

The requirement for both thermal and non-thermal electrons 
does not simply rest on the (rather weak) detection of the source 
at the highest energies in OSSE (3.7(7 and 1.5cr for channels 350- 
600 keV and 600-900 keV, respectively). Instead it is driven by 
the overall spectral shape between ~ 100 and ~ 300 keV where 
the OSSE data have excellent statistics. Even the spectral fits with 
the OSSE data above 350 keV excluded strongly rule out purely 
thermal or purely power-law Comptonization models, as these both 
produce a much steeper high-energy cutoff than observed. 
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Figure 5. PCA, HEXTE and OSSE unfolded data, the best-fitting models and residuals, (a) Comptonization on Maxwellian electrons (TH). (6) Comptonization 
on power-law electrons (PL). Contribution from the multicolour blackbody disc (DISC) and Compton reflection (REF) are also shown. 
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Figure 6. PCA, HEXTE and OSSE unfolded data, the best-fitting models and residuals, (a) with two Comptonizing regions: with power-law (PL) and thermal 
(TH) electrons. (6) Comptonization on hybrid (HYB) electron distribution, (c) hybrid electrons and additional thermal coinponent. Contribution froin the 
multicolour blackbody disc (DISC) and Compton reflection (REF) are also shown. 



5.3.2 Hybrid electrons 

Alternatively, the complex continuum could be produced by a sin- 
gle Comptonizing plasma, where non-thermal electrons partially 
thermalize through Coulomb collisions, creating a hybrid distribu- 
tion: Maxwellian at lower energies and power-law at higher ener- 
gies. We use EQPAIR to compute this distribution self-consistently, 
together with its reflection/line emission and seed photons from the 
disc (model HYB). 

This hybrid model gives a somewhat worse fit ix^ /v = 
100.6/94, see Tableland Fig. |6|)) than the two Comptonizing 
regions model PLTH {x^ jv — 63.4/92), although it is still sta- 
tistically an adequate description of the data. The reason the fit is 
different is that Comptonization from two separate electron distri- 



butions, where one is thermal and the other non-thermal, gives a 
subtly different spectrum to that of a single hybrid electron distri- 
bution. The sum of a thermal and power-law distribution has the 
power-law electrons also contributing to the spectrum at low en- 
ergies, whereas in the hybrid electron distribution there are only 
the thermal electrons at low energies. While in principle this could 
give us a way to distinguish between the two models, in practice 
the situation is much less clear due to uncertainties in the reflected 
spectrum (see below). 

Error ranges are hard to obtain as many of the parameters are 
correlated, so give multiple local minima in x^ . In particular, the 
soft compactness, ^s, is anti-correlated with the fraction of power in 
the non-thermal injection, fnth/^h- This is why in the model HYB 
we fixed = 10. The soft compactness determines the strength of 
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Table 3. Fitting results of the model PLTH, consisting of two Comptoniza- 
tions: on thermal (temperature Ttjj and optical depth T^h) and power-law 
(index p and optical depth Tpi) electrons. Values in parentheses refer to the 
parameters fixed during the fit. The unfolded data and model of this fit are 
shown in Fig.lSU. 
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Compton cooling relative to Coulomb cooling, so a larger fraction 
of the non-thermal electrons thermalize for small l^. The effects 
of this can be offset by putting less power into the thermal heating 
of the electrons i.e. by increasing ^nth/^h- A model HYBN (see 
TableEJ with all the hard power in the non-thermal injection (fixed 
^nth/^h = 1, but is left free) is statistically acceptable (Ax ~ 
+ 1.7 compared to model HYB). This model has is ~ 0.8, which is 
a firm lower limit on the soft compactness. 

An upper limit on the soft compactness can also be found from 
the spectrum. Since the hard-to-soft compactness ratio is very well 
determined by the spectral shape, £h/is ~ 1, large Is requires large 
^h- But there is an upper limit on the hard compactness, since large 
£h gives rise to a strong annihilation line which is not present in the 
data. Thus, we can roughly constrain the soft compactness, I < is 
< 20 and the total compactness, 2 < i < 40. 

The contribution of electron-positron pairs to the Comptoniz- 
ing plasma is negligible, as the best-fitting model yields Tpair ~ 
lO"'^. The plasma would be pair-dominated (rpair > Tion) only for 
extremely high compactness, i > 3000, which is ruled-out by the 
data (x^/f ~ 128/94). 

The Compton reflection requires a substantial amount of rel- 
ativistic smearing. We fix the inner disc radius at _Rrcf ~ WRg 
during the fit, but values higher that ~ 15 Rg are ruled-out. The 
best-fitting model give unphysical rest-frame iron line energy of 
5.6 keV. This probably indicates that the Compton reflection model 
we use here is inadequate, as already inferred from the ASCA fits 
(SectionlsH. 

These fits can be dramatically improved jv — 63.4/92, 
model HYBTH) by including an additional, optically thick, thermal 
Comptonization component (for which we use THCOMP code). The 
residuals in Figs.|6j) and|6l: show that the improvement affects the 
whole spectral range, though in particular the features below ~ 10 
keV are substantially smoothed. The model is detailed in Table|4] 
but again the complex correlations between parameters mean that 
again the errors are only representative as we fixed ^nth/^h = 0.59 
(from model HYB) and E^cf = SOTJg. 

The hot plasma properties are insensitive to the presence of the 
additional thermal Compton component. Both HYB and HYBTH 
models predict £h/^s ~ 1 and a significant fraction of the hard 
power from non-thermal electrons (^nth/^h ^ 0.5). They both rule 
out pair-dominated plasma. The main difference is in the optical 



Table 4. Fitting results with a hybiid thermal/non-thermal model. HYB is 
the hybiid model with fixed = 10 and free ^nth- HYBN is the hybrid 
model with free is and fixed £ntij = 1 (no thermal heating). HYBTH stands 
for the hybrid model with additional component with thermal electrons. 
Values in parentheses refer to the parameters fixed during the fit. The equi- 
librium temperature, Te , was calculated self-consistently and was not a fit 
parameter. 
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depth of the Comptonizing medium, which decreased from ~2.8 to 
~1.8 when THCOMP was added. 

The extrapolated unabsorbed bolometric flux of XTE J 1550- 
564 from both models is 1.1 x 10"'^ erg cm~^ s~^. At a distance of 
5.3 kpc, for a lOM© black hole, this corresponds to about 0.3LEdd- 



6 HIGH AND VERY HIGH X-RAY SPECTRAL STATES 

As we have outlined in the introduction, there are at least two dis- 
tinct soft spectral states of black hole binaries: high and very high. 
We attribute our observation of XTE J1550-564 to the very high 
state. How does it compare with other sources, in both soft states? 
Fig. shows a comparison of the high- and very high-state model 
spectra. The high-state spectrum is from Cyg X-1 (Gierlihski et 
al. 1999) while the very high-state one is from XTE J1550-564 
(this paper). The same model is used for both, with a hybrid elec- 
tron distribution calculated using the EQPAIR code (Coppi 1999). 

The main difference between the high and very high states 
seems to be the amount of Comptonization of the disc spectrum. 
In the high state many of the disc photons escape without being 
strongly Comptonized, therefore the Comptonizing plasma either 
covers only a small fraction of the disc and/or has low optical depth. 
In the very high state most of the disc emission is scattered to sub- 
stantially higher energies, so the Comptonizing medium must have 
a substantial optical depth and cover most of the disc. 

However, we caution that the amount of Comptonization is 
not always straightforward to derive from the spectra. For example, 
Zdziarski et al. (2001) fit two soft state RXTE/OSSE spectra from 
GRS 1915+105. The 'ultrasoft' high-state spectrum (VP 813) has 
a strong, dominating soft component and a weak high-energy tail. 
It is similar to the Cyg X-1 spectrum from Fig.0 but the tempera- 
ture of the soft component and its ratio of power to the hard tail is 
larger. The second observation analyzed by Zdziarski et al. (2001) 
is a typical very high-state spectrum (VP 619), with comparable 
distribution of power between the disc and the tail, similar to the 
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Figure 7. Unabsorbed model components for Cyg X-1 in the high state 
(from Gierlinski et al. 1997) and XTE J1550-564 in the very high state (this 
paper, model HYB). The approximate values of some of the best-fitting EQ- 
PAIR parameters are given. To convert fluxes into luminosities we assumed 
distances of 2 kpc for Cyg X-1 and 5.3 kpc for XTE J1550-564. For Ed- 
dington luminosities we assumed IOMq black hole mass for both sources. 



XTE J1550-564 spectrum analysed here (see fig. 2 in Zdziarski et 
al. 2001). However, the detailed spectral modelling with EQPAIR 
suggests a rather different interpretation of this high state to that 
outlined above. What looks like a disc in the high-state spectrum 
of GRS 1915-1-105 is fit by Comptonization of the disc in optically 
thick, low-temperature plasma, very different to the high state of 
Cyg X-1. However, it seems more likely that this simply reflects 
the fact that the disc spectral models do not adequately describe the 
observed disc emission. The dominating soft component in GRS 
1915-1-105 probably comes from the disc itself, while an optically 
thin, hot plasma is responsible for the weak hard tail. Thus we con- 
clude that there is probably an increase in the covering fraction 
and/or optical depth of the Comptonizing plasma which drives the 
source evolution from the high to very high states. 

Another difference between the high and very high states is a 
change in shape of the non-thermal tail. The power-law index of 
the injected electron distribution is steeper in the very high state 
than in the high state, as is also seen in GRS 1915-1-105 (Zdziarski 
et al. 2001). Thus it seems that there are also changes in the non- 
thermal electron acceleration associated with the source moving 
from high to very high state. 



7 DISCUSSION AND CONCLUSIONS 

During its 1998/1999 outburst the black hole candidate XTE 
J1550-564 showed a variety of X/7-ray spectral states. In this pa- 
per we have analyzed its broadband spectrum during the very high 
state. This spectral state is characterised by the presence of the 
soft, high-energy tail extending without any apparent break above 
~ 300 keV. Such a tail differs from the low/hard state tail both 
in having a significantly softer spectral index and in the lack of 
high-energy cutoff (Grove et al. 1998). Phenomenologically, the 
OSSE spectra of these tails can be described as a power law with 
spectral index F ~ 2.4-3.1 (Grove et al. 1998). In this observation 
r = 3.09 ± 0.04. We have shown that the broadband X/7-ray spec- 
trum of XTE J 1550-564 can be well fit by a model in which seed 
photons from the disc are Compton up-scattered by both thermal 
and non-thermal electrons. Comptonization on purely thermal or 
purely power-law electrons can be ruled out. 

The thermal and non-thermal electron distributions may be 
distinct, or they may instead form a single hybrid thermal/non- 
thermal distribution. This is reflected in the modelling of the XTE 
11550-564 spectrum, which we have fitted by separate thermal plus 
non-thermal (model PLTH) and hybrid (model HYB) Comptoniza- 
tions. However, even if there is a spatially separate non-thermal 
electron population, these electrons will interact via Coulomb col- 
lisions, creating a thermal low-energy extension to the power-law 
distribution. Thus, we expect a hybrid distribution of Comptonizing 
electrons even from a completely non-thermal electron acceleration 
mechanism. Therefore, the PLTH model seems to be less physically 
motivated, and a single hybrid plasma (HYB) would be our model 
of choice as the simplest, physical description of the data. How- 
ever, spectral fitting gives the opposite result, with the less physical 
model (PLTH) being strongly statistically preferred. 

Since we expect that the non-thermal electrons form a hybrid 
plasma, the next step in complexity would be two separate elec- 
tron distributions, with one thermal and one hybrid (rather than one 
thermal and one power law). This model (HYBTH) provides an ex- 
cellent fit to the data, statistically similar to PLTH but without using 
an unphysical power law electron distribution. However, closer ex- 
amination of the Maxwellian component in HYBTH (component 
denoted as TH in Fig.|6l;) reveals its suspicious similarity in spec- 
tral shape to the Compton reflected continuum. Perhaps it is not due 
to emission from some thermal, optically thick plasma at all. Per- 
haps it simply accounts for inaccuracies of the reflection model. In 
particular, the reflection model used here overestimates the depth of 
the iron edge for ionized material as it does not include Compton 
upscattering (Ross, Fabian & Young 1999). While the single hy- 
brid plasma (HYB) is statistically ruled out when compared to the 
two-component solution (HYBTH), better reflection models are re- 
quired before this can be a robust conclusion. 

Whatever is the physical reality behind this additional ther- 
mal Compton component, the properties of the hybrid plasma are 
established quite robustly and do not change much between HYB 
and HYBTH. The ratio of soft to hard compactness is very well 
defined by the shape of the spectrum (£h/is ~ 1), while the frac- 
tion of the power in non-thermal electrons /^h > 0.5. Once 
again, the presence of non-thermal electrons is required to form the 
high energy tail observed beyond 300 keV. But where do they come 
from? 

One obvious mechanism for such high-energy emission is 
from a jet, as in the blazars and radio-loud quasars (see e.g. Ghis- 
ellini et al. 1998). A jet is present in XTE J1550-564 indeed at this 
time, although its radio emission is declining rapidly throughout 
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the OSSE observation (Wu et al. 2002), so its link with the X/7-ray 
emission is doubtful. Moreover, the fact that a similar high-energy 
tail is seen in the high state as well as in the very high state immedi- 
ately rules out a direct association with the jet as the radio emission 
is strongly quenched in the high state (Corbel et al. 2000). 

An alternative mechanism of bulk motion of in-falling mate- 
rial near the black hole horizon was proposed by Chakrabarty & 
Titarchuk (1995), but the free fall results in fairly low electron en- 
ergies, with Lorentz factors 7 < 1.4. These electrons are unable 
to up-scatter photons past ~ 100-300 keV, resulting in a cutoff in 
the power-law spectrum (Laurent & Titarchuk 1999; Zdziarski et 
al. 2001). Since such a cutoff is not seen in the XTE J1550-564 
data, we find the bulk motion Comptonization to be an unlikely 
explanation of the observed emission. 

This leaves magnetic reconnection above the disc as the only 
known possible source of hard X-ray dissipation at high mass ac- 
cretion rates. The accretion disc viscosity is now known to be linked 
to a magnetic dynamo, so magnetic buoyancy can lead to some 
fraction of the reconnection events occurring in a corona rather than 
deep inside the disc (see e.g. Hawley & Balbus 2002), accelerating 
electrons to high energies (di Matteo 1998). We speculate that the 
high and very high state spectra are continuous rather than two sep- 
arate spectral states, linked by an increase in the covering fraction 
of this magnetic corona over the disc. 

Concluding our paper we summarise our main results: 

• The broad-band X/7-ray spectrum of XTE J 1550-564 in the 
very high state can be well described by inverse Compton emission 
from hybrid, thermal/non-thermal plasma. 

• Comptonization on purely thermal or purely power-law elec- 
trons cannot explain the observed spectrum. 

• The high-energy power-law tail observed here is similar in 
other black hole binaries in the high and very high states, show- 
ing that it is a generic property of high mass accretion rate black 
holes. 
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